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ABSTRACT 


Multisensor  track  and  data  fusion  require  accurate  sensor  registration  and  positional  bias  removal. 
This  report  presents  the  problem  of  bias  removal  in  the  context  of  midcourse  interceptor  engagements. 
For  this  application  two  sensors  are  employed:  a  fire-control  radar  and  an  interceptor  with  an  optical 
seeker.  The  fire-control  radar  tracks  objects,  collects  data  on  the  objects,  performs  discrimination,  and 
hands  over  a  target  object  map  to  die  interceptor.  Positional  bias  is  generally  evident  in  the  target  object 
map  and  must  be  removed  in  order  to  fuse  the  radar  and  seeker  data  accurately.  Equations  and  algorithms 
are  derived  for  calculating  the  maximum-likelihood  bias  estimate  for  three  scenarios,  given  three 
assumptions  about  the  measurement  errors  on  the  radar-  and  seeker-tracked  objects.  First,  bias  is  shown 
to  be  an  insignificant  contributor  to  the  assignment  of  objects  when  both  sensors  track  the  same  objects 
and  the  measurement  errors  are  equal  and  uncorrelated  for  all  objects.  Second,  equations  for  direct 
calculation  of  the  optimal  bias  are  derived  for  cases  where  one  sensor  tracks  a  subset  of  the  objects 
tracked  by  the  other  sensor.  Third,  algorithms  are  discussed  for  the  most  complicated  scenarios,  those 
where  dissimilar  objects  are  tracked  by  the  two  sensors. 
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1.  INTRODUCTION 


In  most  ballistic  missile  defense  scenarios,  the  reentry  vehicle  (which  carries  the  warhead)  will  be 
accompanied  by  associated  objects  such  as  tanks  and  deployment  debris.  The  interceptor  (which  has  an 
optical  seeker)  detects  some  of  these  objects  on  its  focal  plane  and  engages  the  most  threatening.  Several 
sources  of  information  may  be  employed  to  correctly  identify  the  reentry  vehicle  as  the  most  threatening 
object  in  the  complex.  In  principle,  discrimination  may  be  accomplished  by  (1)  the  fire-control  radar, 
which  then  hands  over  information  to  the  interceptor;  (2)  the  interceptor  seeker  autonomously,  using  its 
own  measurements;  or  (3)  the  radar  and  seeker  together,  by  combining  data  for  fused  discrimination. 
Regardless  of  the  discrimination  architecture  selected,  handover  of  radar  data  to  the  interceptor  is 
meaningful  only  if  the  objects  seen  by  the  radar  and  those  seen  by  the  seeker  are  correctly  associated  with 
one  another.  This  report  considers  the  handover  of  radar  data  necessary  to  the  interceptor  and  proposes 
robust  algorithms  for  the  association  of  objects  between  the  radar  and  the  seeker.  Figure  1  is  a  cartoon  of 
the  detection  and  handover  process. 


Fire-Control  Radar 


Figure  1.  Handover  of  tracked  objects  from  the  radar  to  the  interceptor. 
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To  ensure  a  successful  handover,  the  information  contained  in  the  radar  uplink  message  must  be 
sufficient  to  unambiguously  associate  the  reentry  vehicle.  A  number  of  signature  and  metric  features  are 
used  for  this  association,  including  angle,  angle  rate,  target  size,  and  scintillation.  The  focus  of  this  report 
is  on  one  of  the  most  commonly  used  correlates,  focal-plane  position.  The  positions  of  the  radar  objects 
are  mapped  to  the  seeker  focal  plane  as  azimuth  and  elevation  positions  and  then  matched  with  the 
positions  of  the  seeker  objects.  This  process  is  referred  to  as  a  “target-object  map”  (TOM).  The  TOM 
may  be  modeled  as  a  linear-programming  problem  and  solved  using  a  minimum-weight  matching 
algorithm  such  as  Munkres  [1,2]  or  JVC  [3,4],  The  linear-programming  problem  searches  for  a  set  of 
variables  that  minimize  a  linear  cost  function  by  taking  into  account  linear  equality  and/or  inequality 
constraints  [5], 

The  generic  assignment  problem  is  shown  in  Figure  2,  where  the  objective  is  to  assign  an  object 
from  one  source  to  an  object  from  another.  The  assignment  problem  is  a  special  subset  of  linear- 
programming  problems  because  it  has  the  constraint  that  all  assignments  are  mutually  exclusive  and 
collectively  exhaustive.  An  example  is  the  assignment  of  a  person  to  a  task,  where  each  person  would  be 
assigned  once  and  only  once  to  a  task  and  each  task  would  be  assigned  once  and  only  once  to  a  person. 
The  constraints  shown  in  Figure  2  ensure  that  each  column  and  each  row  of  the  decision  matrix  D  [where 
D(i,j)  =  djj\  is  assigned  once  and  only  once  and  that  the  decision  variables  are  binary. 


Objective:  minimize 

*  j 

subject  to:  ^dl}  =  l  V  j 

i 

I4  =  1  Vi 

j 

d;j  =  [0  or  1]  V  i,j 

where  dtJ  represents  the  decision  to  assign  (1)  or  not  assign  (0) 

object  i  from  one  source  to  object  j  from  another  source,  and 
is  the  cost  of  the  decision  d;j. 


Figure  2.  Generic  assignment  problem. 
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In  the  TOM  problem  the  costs  cVj  are  derived  from  a  model  of  the  relative  metric  errors  between  the 
two  sensors.  A  few  types  of  errors  can  occur  during  the  uplink  process  that  add  a  translational  bias  to  the 
uplinked  focal-plane  positions  of  the  radar  objects  and  to  the  positions  of  the  seeker  objects.  Random, 
zero-mean,  white  radar  track  and  seeker  measurement  errors  further  complicate  the  association  problem. 
Registration  bias  is  illustrated  on  the  left-hand  side  of  Figure  3;  random  errors  that  result  from 
measurement  noise  are  shown  on  the  right-hand  side.  Measurement  noise  errors  cannot  be  removed; 
however,  bias  can  be  removed  by  shifting  the  uplinked  radar  objects  to  minimize  the  distance  between  the 
radar  and  seeker  objects.  This  optimal  shift  will  be  the  maximum-likelihood  bias  estimate  and  will  be 
embedded  in  the  costs  ci}  shown  in  Figure  2. 


Figure  3.  Bias  error ;  random  error. 


The  objective  of  this  report  is  to  propose  algorithms  for  calculating  optimal  bias.  The  problem  is 
divided  into  three  possible  scenarios.  The  first,  the  mx  n  scenario,  is  defined  with  m  objects  in  common 
between  the  two  sensors.  The  second  scenario,  m  x  n,  is  defined  such  that  m  objects  are  in  common 
between  the  two  sensors  and  one  sensor  has  n  -  m  additional  objects.  The  third  scenario,  the  mismatch, 
assumes  that  at  least  one  object  is  common  and  at  least  one  object  is  different  between  the  two  sensors. 
Venn  diagrams  of  each  scenario  are  shown  in  Figure  4.  The  object  set  for  one  sensor  is  depicted  with 
horizontal  stripes  and  the  set  for  the  other  sensor  has  vertical  stripes,  such  that  the  set  of  objects  in 
common  is  crosshatched. 

The  rest  of  this  report  is  organized  as  follows:  Section  2  defines  the  assignment  problem  and 
develops  the  cost-function  model.  Section  3  develops  four  proposed  algorithms  for  calculating  bias  for 
the  mx  m  and  m  x  n  scenarios.  Finally,  Section  4  offers  a  summary  and  identifies  issues  for  future  work. 
As  the  algorithms  and  discussions  in  this  report  are  applicable  to  problems  beyond  the  radar-to-seeker 
TOM,  the  remainder  of  the  material  will  address  the  problem  generically.  The  reader  is  encouraged  to 
draw  parallels  with  other  two-sensor  tracking  correlation  problems  where  registration  biases  play  an 
important  role. 
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mxm  mxn  Mismatch 


Figure  4.  Venn  diagrams  of  scenarios. 
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2.  STATEMENT  OF  PROBLEM 


The  objective  of  the  track-to-track  assignment  problem  is  to  find  the  optimal  minimum  error 
assignment  of  objects  viewed  by  two  sensors.  This  problem  is  complicated  by  the  presence  of  registration 
biases  that  shift  the  positions  of  all  objects  from  truth.  The  total  bias  may  have  an  effect  on  the 
assignment;  therefore,  bias  must  be  accounted  for  to  ensure  the  accuracy  of  the  track-to-track  assignment. 
The  generic  assignment  problem  shown  in  Figure  2  will  now  be  tailored  using  the  particular  features  that 
distinguish  the  track-to-track  assignment  problem,  such  as  bias  and  the  number  of  common  objects. 

Bias  is  defined  as  the  shifted  difference  in  position  between  the  objects  of  the  two  sensors.  Let 
be  the  azimuth  and  elevation  positions,  respectively,  of  the  zth  object  for  sensor  one;  let  be  the 
azimuth  and  elevation  positions,  respectively,  of  the  7th  object  for  sensor  two.  If  there  are  no  random 
measurement  errors,  bias  [bu  bv),  affects  all  targets  equally  so  that  uf  =  x,  -  bu  and  v;  =  yt  -  bv  when  4  =1. 
Given  the  presence  of  random  errors,  the  optimal  assignment  may  be  calculated  by  weighting  the 
positional  differences.  When  the  errors  are  Gaussian,  the  weights  are  the  measurement  error  covariance 
matrices  for  the  zth  and/th  objects  of  sensors  one  and  two,  P,  and  2},  respectively.  The  cost  function  from 
Figure  2,  rewritten  to  include  the  assumed  error  model,  closely  resembles  a  linear  least-squares  cost  as 
shown  in  Figure  5.  The  cost  function  in  Figure  5  minimizes  the  log  likelihood,  as  derived  in  the 
Appendix. 


/I 

f  X-b„  -U:  1 

T 

/  \  — 1 

X;  -bu~U;~ 

c„  =  m(|/> 

)+HM)+ 

1 

1 

* 

1 

(^j 

l  U  J 

_y>-bv-Vj- 

where  [x,-  y,]  is  the  position  of  the  zth  object  for  sensor  one, 

[uj  Vj]  is  the  position  of  the 7th  object  for  sensor  two, 

[bu  bv]  is  the  bias, 

Pt  is  the  covariance  matrix  of  the  zth  object  for  sensor  one,  and 
Zj  is  the  covariance  matrix  of  the 7'th  object  for  sensor  two. 


Figure  5.  Weighted-distance  cost  Junction. 
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Given  the  presence  of  bias,  the  problem  is  nonlinear  programming  with  respect  to  bu.  bv,  and  (dl}) 
instead  of  a  linear  problem  of  (d,f)  alone  (as  in  Figure  2).  Ideally,  an  attempt  would  be  made  to  solve  the 
problem  presented  in  Figure  5  directly,  solving  for  bias  and  decision  variables  at  the  same  time;  however, 
it  is  difficult  and  computationally  expensive  to  solve  a  nonlinear  binaiy  problem  accurately.  Alternative 
algorithms,  therefore,  have  been  developed  that  calculate  bias  before  solving  for  decision  variables;  these 
alternative  algorithms  are  presented  in  this  report. 

The  constraints  (see  Figure  6)  also  change  from  those  shown  in  Figure  2  to  allow  for  an  unequal 
number  of  objects  between  the  two  sensors.  (Note  that  if  both  sensors  track  the  same  objects,  the 
constraints  in  Figure  6  are  the  same  as  those  in  Figure  2.)  Without  loss  of  generality,  it  is  assumed  in  this 
report  that  m  <  n;  it  is  also  assumed  that  each  sensor  contributes  all  the  objects  that  it  is  able  to  track  to 
the  track-to-track  assignment  problem. 


Assumptions;  ( 1 )  sensor  one  tracks  m  objects 

(2)  sensor  two  tracks  n  objects,  where 

(a)  objects  [1 :  m]  are  the  same  as  the  m  objects  for  sensor  one 

(b)  objects  [m  + 1 :  n]  are  extra  objects 

Constraints: 

X4  =  l  V;'€[l:m] 

i 

5X‘=0  V  j  e[m  +  l:n] 

i 

2X  =  1  V  i  e  [1 :  m] 


Figure  6.  Constraints  adjusted for  an  unequal  number  of  objects. 


Although  the  most  general  form  of  the  bias  calculation  problem  can  be  complex  and 
computationally  expensive  to  solve,  assumptions  can  be  made  to  constrain  the  problem  and  improve  its 
solvability.  These  special  cases  provide  valuable  insight  into  the  track-to-track  assignment  problem.  For 
each  of  the  three  scenarios,  m  x  m,  m  x  «,  and  mismatch,  one  of  three  assumptions  may  be  made  about  the 
error  covariances  of  the  objects  for  a  given  sensor.  As  mapped  in  Table  1,  five  different  bias  calculation 
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algorithms  have  been  developed  to  handle  these  three  scenarios,  as  well  as  the  three  assumptions.  First,  it 
may  be  assumed  that  the  covariance  matrices  for  a  given  pair  of  objects  sum  either  to  the  identity  matrix, 
P’i  +  'Lj  =  I,  or  to  a  diagonal  matrix,  P  +  £;  =c/,  where  c  is  any  vector  [<r*  ].  (Note  that  cost 

function  is  reduced  to  a  Euclidean-distance  function  when  P  +  £  .  =  /  .)  Second,  it  may  be  assumed  that 
the  covariance  matrices  for  each  sensor  are  uncorrelated  and  identical  for  all  objects  such  that 
P(  +  Xj i  =  W,  where  W  denotes  a  constant  matrix  that  is  independent  of  /  or  j.  Third,  if  no  assumptions 
are  made,  the  covariance  matrices  are  modeled  as  unique  for  each  object  such  that  P  +  X ;  =  Wip  where 
^denotes  a  matrix  that  is  dependent  on  i  and  j.  The  remaining  sections  of  this  report  present  algorithms 
for  calculating  bias  for  these  three  scenarios  and  these  three  assumptions.  Although  real-time  algorithms 
will  not  be  able  to  use  scenario  typing,  stepping  through  these  algorithms  for  special  cases  can  assist  in 
the  formulation,  albeit  suboptimal,  of  bias  removal  algorithms  for  the  general  problem. 


TABLE  1 

Map  of  Assumptions,  Scenarios,  and  Proposed  Algorithms 


Covariance 

Assumption 

Scenario 

Assumption 

Case 

Bias  Algorithm* 

+  + 

II  II 

a,  ^ 

mx  m 

1 

Insignificant 

mxn 

3 

Loop  through,  Equation  (10) 

mismatch 

5 

Loop  through,  Equation  (7) 

mx  m 

2 

Direct  estimation  Equation  (8) 

Pi+lj=  w 

mx  n 

3 

Loop  through,  Equation  (10) 

mismatch 

5 

Loop  through,  Equation  (7) 

P/-+X/  -  Wjj 

mxm 

4 

Loop  through,  Equation  (6) 

mx  n 

4 

Loop  through,  Equation  (6) 

mismatch 

5 

Loop  through,  Equation  (6) 

‘Equations  are  shown  in  Section  3. 


3.  BIAS  CALCULATION  ALGORITHMS 


3.1  CASE  1 :  BIAS  AS  AN  INSIGNIFICANT  CONTRIBUTOR 

This  first  case  is  presented  in  order  to  understand  the  effect  of  bias  on  the  optimal  assignment  for 
the  m  x  m  scenario  when  P.  +  X , .  =  cl.  It  will  be  shown  that  bias  may  be  analytically  factored  from  cost 
function,  thereby  making  bias  an  insignificant  contributor  to  the  optimal  decision.  To  understand  how 
bias  affects  cost  function,  the  objective  function  from  Figure  5  is  first  manipulated.  (The  problem  is 
worked  through  using  only  scalar  positions  in  order  to  preserve  the  simplicity  of  the  mathematics.) 


n  ;  ;  ' 


Cj  i  j 


Cl 

J_ 

Cl 


(1) 


The  following  constraints  for  an  m  x  m  problem  are  assumed: 


(2) 


Inserting  the  constraints  from  Equation  (2)  into  (1)  gives  the  following  objective  function: 


z- 


u)  +  him - 2 bu 'Lxi  - 2 bu 


\ 


(3) 


The  objective  in  the  assignment  problem  is  to  select  values  for  decision  variables  such  that  cost 
function  is  minimized.  If  a  term  in  the  cost  function  does  not  depend  on  the  decision  variables,  its 
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contribution  to  the  cost  of  the  decision  is  that  of  a  constant:  that  is,  it  does  not  encourage  or  discourage  a 
given  selection  of  the  decision  variables.  Employing  this  observation  simplifies  the  objective  function  in 
Equation  (3)  to  z\  where  z  -  ^  +  k2z,  and  k,  and  k2  are  constants,  as  follows: 


—  Zxf  +  Zu2  +  b2um  -  2bu  Z -  2 bu  Zuj 
C\  V  1  J  i  j 

c 

z  -  zr(z  -  constant)  =  -ZZx^d 

1  i  j 


=  constant 

(4) 


Bias  does  not  appear  in  the  new  objective  function  in  Equation  (4)  and  therefore  does  not  influence 
the  solution.  The  optimal  assignment  of  objects  is  independent  of  bias  in  the  m  x  m  scenario,  where  both 
sensors  track  the  same  objects  and  the  measurements  errors  are  uncorrelated  i^P  +  Zj  =  clj.  The 
algorithm  for  this  first  case  thus  requires  no  special  treatment  of  bias  prior  to  applying  one'  of  the  classic 
assignment  algorithm  solvers  such  as  Munkres  or  JVC. 

3.2  GENERAL  SOLUTION  OF  BIAS 

Although  no  bias  removal  is  necessary  for  the  m  x  m  scenario  when  P  +  =cl,  it  is  necessary  in 
all  other  cases  that  are  to  be  considered.  A  general  equation  for  calculating  bias  is  derived  and  discussed 
in  this  section  and  is  referenced  in  the  remaining  sections.  To  solve  for  the  bias  shown  in  Figure  5,  the 
objective  function  is  manipulated  as  follows: 


z=ZZ\ 

*  j 

dz 


InfWMM 


+ 


-  uj  -  K 


IT 


y<  -  Vj  -  K 

a 


£A_a_ 

db 


xi  ~  uj  ~  b 

>\  -  v,  -  b 


Xt  -  Uj  -  bu 


V*,r 


n  \ 


xi  ~  uj  -  bu 

[y.-yj-Kp 


(5) 


a?  m  "  /  .  X;  —  U:  -b, 


V,  . 


The  equation  for  bias  may  be  formed  by  setting  dz  /  db  =  0  and  then  solving  for  b : 
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(6) 


■Uj 


m  n  /  x_i 

ds-SX/>+Z; 

■  ; 


d;; 


xi-uj 


,y<  ~  vj . 


1  A 
d- 

y 


The  general  solution  for  the  maximiun-likelihood  bias  estimate  is  shown  in  Equation  (6).  Bias  cannot  be 
calculated  directly  using  Equation  (6),  as  the  unknown  decision  variables  (di})  are  embedded.  The  error 
covariance  assumptions  mapped  out  in  Table  1  are  employed  in  the  remaining  sections  of  this  report  to 
manipulate  Equation  (6)  to  be  a  more-applicable  estimator. 

3.3  CASE  2:  DIRECT  CALCULATION  OF  BIAS 

Case  2  is  defined  by  the  assumption  that  Pi  +  'Lj=cI  or  P.+  £.  =  W  for  the  m  x  m  scenario. 
Although  it  has  already  been  proven  that  bias  is  an  insignificant  contributor  when  P  +  =  cl  for  this 

scenario  it  may  be  desirable  to  calculate  it.  The  general  solution  for  bias,  shown  in  Equation  (6),  can  be 
manipulated  as  follows  for  case  2: 
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Equation  (7)  may  be  manipulated  by  employing  the  constraints  identified  in  Equation  (2),  such  that 
the  equations  for  the  optimal  bias  are: 
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I*,  Ly,  -Sv; 

ft.  =  J - J—  and  h  =  -i - i— 


m 


m 


Bias  is  independent  of  the  decision  variables  and  straightforward  to  compute  using  Equation  (8); 
therefore,  bias  may  be  directly  calculated  when  Pt  +  X,  =  c/  or  P:  +  I  =W  for  the  m  x  m  scenario.  This 
solution  is  the  well-known  centroid  method  for  bias  removal.  Bias  is  the  difference  in  the  positional 
centroids  for  the  objects  of  each  sensor.  This  method  is  optimal  if  the  given  assumptions  hold  true; 
however,  this  method  may  be  significantly  suboptimal  if  these  conditions  do  not  hold  true,  especially  the 
assumption  that  each  sensor  tracks  the  same  m  objects. 

3.4  CASE  3:  PROCESS  POSSIBLE  COMBINATIONS 

The  problem  increases  in  difficulty  when  the  m  x  n  scenario  is  considered.  In  the  m  x  n  scenario, 
one  sensor  tracks  m  objects  and  the  other  tracks  n  objects,  where  the  m  objects  are  a  subset  of  the  n 
objects.  As  there  are  an  unequal  number  of  objects,  bias  is  a  significant  contributor  to  the  overall 
assignment,  even  when  P  +  .  =  cl  is  assumed.  The  general  solution  for  bias,  shown  in  Equation  (6), 

can  be  manipulated  as  follows  for  case  3: 


v  = 

A. 


f  m  n 

=  (  lW~'d, 


\  l(  m  n 


V*  j 


IIW 

i  j 


-1 


*,  -Uj 


y,-vj 


\d, 


1  j 


( 

m  n 

X; 

\ 

ww~l 

II 

l1'  ; 

1 

I 

dij 

j 

m 


1Lxitldij-'IJUjt,dij 
‘  j  j  ' 

m  n  n  m 

ly^-lvjld, 


(9) 


The  algorithm  proposed  for  case  3  is  used  to  calculate  the  bias  and  cost  of  the  optimal  assignment 
for  a  given  combination  of  m  objects.  The  idea  behind  this  method  is  that  if  one  sensor  tracks  two  objects 
and  the  other  tracks  three,  three  candidate  matches  can  be  searched  for  the  optimal  match,  as  there  are 
n!/m!(n-m)!  possible  ways  of  choosing  m  out  of  n  objects.  Choosing  objects  1  and  2  means  that  object  3 
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from  the  second  sensor  would  not  be  assigned  to  one  of  the  objects  from  the  first  sensor.  The  constraints 
for  the  m  x  n  case  are  defined  in  Figure  7. 


Let  Sk,  k  =  1:  — — ,  be  a  set  of  integers  {;'</,<•••<  jm  <«)  where  n>m. 

m\(n-m)\  1  J 

For  each  k,  the  following  constraints  are: 

X4/  =1>  ^  ~  0,  and  XX4  =  1 

‘  jest  i  / 


Figure  7.  Constraints  for  the  mXn  scenario. 


The  constraints  shown  in  Figure  7  are  used  in  conjunction  with  Equation  (9)  to  derive  the  following 
equations  for  bias: 


Xx,  -  X  k-  X>',  -  X  v, 

fru=J - and  &v  =  - - J^— 

m  m 


Therefore,  no  direct  equation  exists  for  bias  that  is  independent  of  the  decision  variables  for  the  m  x  n 
scenario  when  Pt  +  X^ ,  =  cl  or  Pt  +  X7  =  W .  The  bias  that  is  calculated  for  a  given  combination  is  used 
as  input  into  the  objective  function,  and  a  track-to-track  assignment  is  calculated  for  that  combination. 
This  procedure  is  performed  for  each  possible  combination,  and  the  assignment  with  the  smallest  cost  is 
chosen  as  the  optimal  assignment. 

3.5  CASE  4:  BIAS  REMOVAL  WHEN  COVARIANCES  ARE  DISSIMILAR 

The  algorithm  for  case  4  addresses  the  bias  calculation  problem  for  the  m  x  m  scenario  and  the 
m  xn  scenario  when  /*  +  Xy  =  WtJ  such  that  WtJ  depends  on  i  and  j.  Direct  equations  (8)  and  (10)  for  bias, 
*  were  derived  from  Equation  (6)  however,  when  P  +  =  W,  the  decision  variables  cannot  be  removed 

from  Equation  (6),  and  there  is  no  explicit  equation  for  calculating  bias  independently  of  decision 
variables.  To  handle  this  case,  another  loop  is  added  to  the  algorithm  such  that  each  possible  assignment 
of  m  objects  is  processed,  and  for  each  assignment  the  bias  is  calculated  using  Equation  (6).  For  each 
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iteration,  the  optimal  track-to-track  assignment  is  determined  using  the  calculated  bias.  The  overall 
optimal  assignment  is  then  the  solution  with  the  smallest  cost. 

As  there  generally  will  be  more  possible  than  feasible  assignments  to  process,  adaptations  to  this 
algorithm  can  significantly  decrease  the  number  of  iterations.  One  technique  is  to  use  the  measurement 
error  covariance  matrix  to  determine  whether  an  assignment  is  feasible.  The  m  x  m  scenario  has  only  a 
few  feasible  assignments  where  all  objects  from  one  sensor  are  assigned  to  the  objects  from  the  second 
sensor  and  are  within  a  given  bound  of  each  other.  If  an  error  bound  on  the  magnitude  of  the  bias  is 
available,  it  can  also  be  used  to  determine  feasible  assignments.  An  assignment  is  assumed  to  be  credible 
if  before  bias  removal  the  objects  from  the  two  sensors  are  within  a  given  threshold  of  each  other  such 
that 


Ti 


uj 


-  V: 


(Pi+ij+PtY 


<;i2 . 


(ii) 


An  example  of  these  bounds  is  given  in  Figure  8,  where  A,  B,  and  C  are  the  objects  from  one  sensor,  and 
1, 2,  3,  and  4  are  the  objects  from  another  sensor. 


Figure  8.  Elimination  of  possible  combinations  using  the  bias  error. 


A  Boolean  combination  matrix  may  be  constructed  that  identifies  all  possible  initial  combination  options. 
[See  Equation  (12)  for  the  example  in  Figure  8.]  Note  that  a  value  of  one  means  the  assignment  is 
feasible;  zero  means  it  is  infeasible. 
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(12) 


c= 


10  0 
0  1  0 
_0  0  1 
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This  combination  matrix  reduces  the  number  of  iterations  from  24  to  2.  The  idea  is  to  pattern  match  for 
objects  that  can  be  reasonably  assigned;  for  example,  there  would  be  no  need  to  try  to  match  objects  A 
and  3  from  Figure  8. 

Another  adaptation  is  to  assume  that  one  specific  object  will  always  be  detected  and  discriminated 
by  sensor  one  and  detected  by  sensor  two,  therefore  assuming  that  the  object  will  definitely  be  assigned 
from  each  sensor-object  set.  Another  assumption  that  improves  the  accuracy  of  the  calculated  bias  and 
decreases  computational  complexity  is  that  each  sensor  detects  and  discriminates  one  object 
unambiguously;  bias  can  easily  be  calculated  as  the  difference  between  the  positions  of  these  two  objects. 
Explanation  and  development  of  further  adaptations  will  be  presented  as  future  work. 

3.6  CASE  5:  COMMENTS  ON  THE  MISMATCH  SCENARIO 

Calculating  bias  and  solving  for  optimal  decision  variables  is  a  complex  task  when  the  two  sensors 
track  different  objects.  An  example  of  the  mismatch  scenario  is  one  sensor  tracking  five  objects  and 
another  sensor  tracking  four,  but  only  three  of  the  objects  are  common.  The  number  and  identity  of 
common  objects  is  unknown  a  priori,  as  well  as  the  correct  value  of  the  bias.  The  track-to-track 
assignment  problem  for  the  mismatch  scenario  is  nontrivial  because  of  this  list  of  unknowns.  Without 
making  assumptions,  however,  this  is  the  real-time  scenario;  therefore,  algorithms  must  be  robust  to 
operate  in  the  mismatch  scenario. 

The  basic  approach  to  solving  the  problem  proposed  here  expands  on  the  algorithm  presented  in 
case  3  for  solving  the  m  x  n  scenario  when  P,  +  X,  =  W .  In  the  m  x  n  scenario,  m  objects  are  chosen 
from  the  n  objects  and  used  to  calculate  a  bias  and  a  track-to-track  assignment.  The  assignment  with  the 
smallest  cost  is  chosen  as  optimal.  A  similar  approach  is  followed  for  this  case,  except  now  the  number 
of  possible  combinations  is  significantly  large.  It  is  no  longer  known  that  there  will  be  m  matched 
objects,  so  their  number  will  be  varied  from  1  to  m,  with  all  possible  combinations  processed. 

In  this  multiple-objective  assignment  problem,  the  aims  are  to  (1)  calculate  bias  and  find  decision 
variables  to  minimize  the  total  distance  of  all  assigned  objects  and  (2)  maximize  the  total  number  of 
assigned  objects.  Two  possible  solution  algorithms  are  discussed.  In  the  first  the  two  goals  are  weighted 
to  create  the  new  objective  function  shown  in  Figure  9.  Unfortunately,  the  weights  are  problem- 
dependent  parameters  chosen  through  experience.  If  the  only  goal  were  to  minimize  the  distance,  the 
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solution  would  always  be  to  assign  only  one  object  such  that  bu  =  x,  -  u},  bv  =  y,  -  y,,  and  z  =  0.  If  the  only 
goal  were  to  maximize  the  number  of  assigned  objects,  m  objects  would  always  be  assigned,  which  does 
not  necessarily  give  the  correct  pattern  match. 

The  second  algorithm  assumes  that  each  cost  zk  is  a  chi-square-distributed  random  variable  with  nk 
degrees  of  freedom  (ndof),  such  that  ndof  =  2  x  #  objects  assigned.  Because  each  solution  may  have 
different  degrees  of  freedom,  each  solution  must  be  converted  to  have  a  common  degree  of  freedom  in 
order  to  compare  costs.  An  algorithm  to  map  the  cost  zk  with  nk  degrees  of  freedom  to  zk  with  one  degree 
of  freedom  is  shown  in  Figure  10.  The  costs  for  each  possible  assignment  are  adjusted  for  the  number  of 
objects  assigned,  and  the  solution  is  simply  chosen  as  the  minimal  zk-  This  algorithm  does  not  highly 
encourage  solutions  where  many  objects  are  assigned.  To  more  strongly  encourage  the  second  objective, 
extra  weighting  may  be  applied,  similar  to  the  first  algorithm. 


Objective:  maximize  v  =  ^  (cxqk  -  c2zk  )dk 

k 

subject  to:  dk  -  [0  or  1]  V  k 

2X= i  v* 

k 

where  qk  is  the  number  of  objects  assigned  for  the  kth  iteration, 
zk  is  the  optimal  cost  of  the  kth  iteration  after  bias  removal, 
q  is  the  weight  of  the  goal  to  maximize  the  number  of  objects  assigned, 
c2  is  the  weight  of  the  goal  to  minimize  the  cost  of  the  assignment,  and 
dk  is  the  decision  to  use  the  solution  of  the  Mi  iteration. 


Figure  9.  Multiple -objective  assignment  problem. 
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Let  x  be  a  %2  distributed  random  variable  with  n  degrees  of  freedom. 
Let  x'  be  a  %2  distributed  random  variable  with  one  degree  of  freedom. 
Let/be  the  cumulative  distribution  function  that  maps  a  x2 
random  variable  with  n  degrees  of  freedom  to  a  cumulative 
probability,  p,  such  that  p  =  f(x,n). 

Let  g  be  a  function  that  is  the  inverse  function  off,  such  that 
x  =  g{p,n )  =  g(f(x,n),n). 

Therefore,  x  is  mapped  to  x'  such  that  x'  =  g(f(x,n),  l). 


Figure  10.  Algorithm  to  map  a  cost  to  one  degree  of  freedom. 


* 


* 


The  proposed  algorithms  for  the  mismatch  scenario  are  not  necessarily  practical  for  a  large  number 
of  objects,  as  the  number  of  iterations  increases  rapidly.  For  example,  if  sensor  one  tracks  three  objects 
and  sensor  two  tracks  three  objects,  there  are  nine  ways  to  assign  one  object,  nine  ways  to  assign  two 
objects,  and  one  way  to  assign  three  objects,  which  totals  19  iterations.  Total  iterations  for  a  problem  with 
m  and  n  tracked  objects  are: 
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*=1  v 


ml 
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k\(m-k)l 


n\ 


k\{n-k)\ 


The  number  of  iterations  discussed  in  this  section  are  for  processing  the  possible  combinations  when 
Pi  +  'Lj=  W.  When  Pi  +  'Lj=  W~  significantly  more  iterations  must  be  processed  in  addition  to  each 
possible  assignment  discussed  in  this  section.  The  adaptations  discussed  in  Section  3.5  are  equally 
applicable  to  decreasing  the  number  of  iterations  for  the  mismatch  scenario. 

3.7  SUMMARY  OF  CASES 

The  algorithms  described  in  cases  2  through  4  are  related  as  shown  in  Figure  1 1.  For  case  2,  the 
m  x  m  scenario  where  P  +  X;  =  cl,  bias  may  be  directly  calculated  using  Equation  (8).  Case  3  is  similar 
to  case  2  except  that  the  m  x  m  scenario  is  assumed,  so  bias  is  calculated  by  iterating  through 
Equation  (10)  for  all  possible  combinations  of  objects.  For  case  4  an  extra  loop  is  added  to  the  flowchart 
in  Figure  11,  as  all  possible  assignments  are  processed  using  Equation  (6).  In  the  mismatch  scenario, 
case  5,  Equation  (7)  is  cycled  through  to  calculate  the  bias  for  P  +  =  W ,  and  Equation  (6)  is  used  for 

P  +  ’Lj  =  Wfj.  As  the  number  of  iterations  grows  rapidly  in  the  mismatch  scenario,  even  for  a  small 
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number  of  objects,  it  is  necessary  to  use  techniques  to  identify  the  feasible  iterations  to  reduce  computer 
runtime. 


Figure  11.  Flowchart  of  algorithms  for  cases. 


So  far  this  report  has  emphasized  calculating  bias  as  an  unknown  parameter;  however,  after  bias  is 
calculated  it  can  no  longer  be  modeled  in  the  cost  function  as  in  Figure  5.  With  a  given  mean  ( ba ,  bv)  and 
an  error  covariance  Pb,  bias  is  modeled  as  a  Gaussian  random  variable.  The  new  cost  function  to  use  in 
solving  for  the  decision  variables,  is  shown  in  Figure  12. 
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c.  =ln(|^l)  +  l„(K|)+ln(|Zj|)+  (p.  +  n  +  ij)-1  *‘-_bruj 

y  i  u\  Vj  Ji  U\  vj 


where  ]bu  £v]  is  the  bias  estimate  and  Pb  is  the  error  on  the  bias  estimate 


Figure  12.  Cost  function  using  estimated  bias. 
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4.  CONCLUSIONS 


Accurately  calculating  and  removing  bias  is  pertinent  to  the  accurate  association  of  objects  between 
two  sensors.  This  report  presents  five  algorithms  for  calculating  bias.  The  algorithms  cover  three 
scenarios  and  three  assumptions  that  are  based  on  the  number  of  similar  objects  and  measurement  errors, 
respectively,  for  the  radar  and  seeker.  In  a  few  simple  cases  bias  was  shown  to  be  negligible  or  easily 
calculated.  It  was  proven  that  for  case  1  bias  is  an  insignificant  contributor  to  the  optimal  assignment  of 
objects.  Two  equations  were  derived  to  calculate  directly  the  optimal  bias  for  cases  2  and  3.  For  cases  4 
and  5  it  was  shown  that  bias  is  more  difficult  to  calculate  than  for  other  special  cases  and  may  require  an 
iterative  algorithm.  Advanced  adaptations  of  the  algorithms  were  briefly  discussed;  they  can  drastically 
help  reduce  the  number  of  required  iterations  in  many  cases.  Robust  algorithms  that  use  these  adaptations 
will  be  studied  further  and  implemented  as  future  work. 


APPENDIX 

DERIVATION  OF  WEIGHTED-DISTANCE  COST  FUNCTION 


The  cost  function  in  Figure  2  is  replaced  in  Figure  5  by  a  standard  weighted-distance  cost  function. 
Its  derivation  is  an  adaptation  of  a  derivation  of  the  weighted-distance  cost  function  that  assumes 
P  +  'Lj  =  W  [6],  The  objective  is  to  find  the  maximum-probability  hypothesis  Hy  that  object  i  from 
«  sensor  one  is  associated  with  object  j  from  sensor  two.  The  maximum-likelihood  estimate  of  the  true 

position  of  the  object  is  £;  xi  is  the  estimate  from  sensor  one;  and  is  the  estimate  from  sensor  two. 
This  derivation  is  for  a  one-dimensional  distance  instead  of  the  standard  two-dimensional  (azimuth  and 
*  elevation)  to  simplify  the  mathematics. 


maximize  pi}  -  P[H -) 

Pa  =  p{xi\z)p{uM)  (A_1) 


Maximizers  of  py  also  maximize  any  monotone  function  of  py.  As  the  natural  logarithm  is  a  monotone 
function,  the  maximization  problem  of  Equation  (A.  1)  may  be  restated  as  follows: 


maximize pi}  <=>  minimize  c{j  =  -ln(Py)  • 


(A-2) 


Assuming  a  Gaussian  distribution  function  for  the  position  of  the  object  in  Equation  (A.  3),  the 
following  simplifications  of  Equations  (A.4)  through  (A.7)  may  be  performed: 


Pij  = 


-z)7  pi~'  (*t  -z)+(“;  V  (“j  "t)j 


(A-3) 
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(A-4) 


(A-5) 


(xi -Z?P~\*i  ~z)  =  (xi~uj)  +1;)"^^  +  I,-)_1(jCi  ~Uj) 

(uj  - zf  1  -l(uj  - z)  =  (xt  - Uj)T(p,  + 1,)'1 1  .(/>  +  X,.)>;  - Uj)  ’ 


(A-6) 


(Xi  -  zf  P^(X.  -  z)  +  (“,  -  z)  l~l(uj  -z)  =  (x;  -  «/(/>  +  Z,.)"1^  -Uj)  .  (A-?) 


Using  Equation  (A. 7),  (A.3)  becomes 


n  = 


(A-8) 


The  weighted-distance  cost  function  is  therefore 


cv 


(A-9) 


for  the  one-dimensional  position  case.  The  two-dimensional  position  equivalent  of  (A.9)  is  shown  in 
Figure  5. 
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Multisensor  track  and  data  fusion  require  accurate  sensor  registration  and  positional  bias  removal.  This  report 
presents  the  problem  of  bias  removal  in  the  context  of  midcourse  interceptor  engagements.  For  this  application  two 
sensors  are  employed:  a  fire-control  radar  and  an  interceptor  with  an  optical  seeker.  The  fire-control  radar  tracks 
objects,  collects  data  on  the  objects,  performs  discrimination,  and  hands  over  a  target  object  map  to  the  interceptor. 
Positional  bias  is  generally  evident  in  the  target  object  map  and  must  be  removed  in  order  to  fuse  the  radar  and  seeker 
data  accurately.  Equations  and  algorithms  are  derived  for  calculating  the  maximum-likelihood  bias  estimate  for  three 
scenarios,  given  three  assumptions  about  the  measurement  errors  on  the  radar-  and  seeker-tracked  objects.  First,  bias 
is  shown  to  be  an  insignificant  contributor  to  the  assignment  of  objects  when  both  sensors  track  the  same  objects  and 
the  measurement  errors  are  equal  and  uncorrelated  for  all  objects.  Second,  equations  for  direct  calculation  of  the 
optimal  bias  are  derived  for  cases  where  one  sensor  tracks  a  subset  of  the  objects  tracked  by  the  other  sensor.  Third, 
algorithms  are  discussed  for  the  most  complicated  scenarios,  those  where  dissimilar  objects  are  tracked  by  the  two 
sensors. 
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